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A developing surface frontal system with an overlying polar front 
jet and subtropical jet stream merging region was present over the con- 
tinental United States during 22 to 24 March, 1966. Possible cross- 
over flow between these two jet streams has been studied in detail by 
isentropic air trajectories at four potential temperature levels from 
310 K to 335 K. 

The results of the trajectory analyses on the four 0 surfaces 
have indicated a cross-over of air parcels from the subtropical jet 
stream to the polar front jet stream. No evidence has been found for a 
reciprocal relationship for air parcels initially associated with the 
polar front jet between 310 K and 335 K isentropic levels. 
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SECTION I 



INTRODUCTION 

Studies by Reiter and Whitney (1965) of the region where a polar 
front jet and a subtropical jet stream are in close proximity to each 
other have produced a provocative hypothesis that there exists a cross- 
over flow between the two respective jets. One corollary is that cer- 
tain air parcels in the polar front jet, rather than follow a normal 
cyclonic trajectory over the continental United States, cross-over and 
proceed along a trajectory path similar to that of air parcels in the 
subtropical jet stream, A second corollary proposes that, rather than 
follow a basically anticyclonic trajectory, certain air parcels in the 
subtropical jet cross-over and then continue along a trajectory path 
similar to that of air parcels in the polar front jet stream flow. 

Regardless of the validity of the hypothesis, the isentropic tra- 
jectory of a representative air parcel in either jet can be computed by 
the method developed by Danielsen (1961). This computational method 
describes with considerable accuracy the movement of air parcels along 
an isentropic surface. However, unlike several other approaches to 
calculating air trajectories which assume that air parcels move quasi- 
horizontally on a surface of constant height or pressure, the Danielsen 
approach gives the necessary third dimension — a measurement of the 
average vertical motion. 

Throughout this study diabatic effects on trajectory motion in 
this non-precipating region are neglected because diabatic contribu- 
tions to vertical motion are normally small compared to adiabatic con- 
tributions . 
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During late March 1966 a strong developing surface low with as- 
sociated frontal systems moved east-northeast across the continental 
United States. Also associated with this surface system were a strong 
polar front and a subtropical jet stream in close proximity to each 
other as illustrated in Figure 1. The period between 22 March 0000 GMT 
to 24 March 0000 GMT was selected as being most favorable for testing 
the hypothesis of Reiter and Whitney (1965). 
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Fig. 1. 23 March 1966, 1200 GMT subtropical and polar front jet stream analysis on the 300 mb 

isobaric surface. Dashed lines are isotachs (units kts). 



SECTION II 



NATURE AND ACQUISITION OF THE DATA 

All the data in the analysis comes from radiosonde and rawinsonde 
checked data found in the U-S, Weather Bureau Northern Hemisphere Data 
Tabulations for 22, 23, and 24 March 1966. 

The period under study contains five observation times, 22 0000 
GMT, 22 1200 GMT, 23 0000 GMT, 23 1200 GMT, and 24 0000 GMT, For each 
observation time, an area of interest is determined by plotting an area 
encompassing not only the surface system but also the region ot the 
merging jets. Table 1 shows the data areas for each observation time. 
Having delineated the regions of interest, the checked data is utilized 
from each upper-air station within these regions. 

The data for each upper-air station consists of all standard pres- 
sure levels in 50 mb intervals from 1000 mb to 200 mb. From 200 mb to 
100 mb the interval is 25 mb. Temperature in degrees Celsius, heights 
in meters, wind direction in degrees, and wind speed in meters per sec- 
ond are tabulated for each standard level. In addition to the standard 
levels, the significant levels are recorded with their tabulated temp- 
eratures and heights 
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Table 1. Data areas used for each observation time 



OBSERVATION 

TIME 


NORTH 


SOUTH 


WEST 


EAST 


22 March 
0000 GMT 


U. S. -Canada 
Border 


U.S, -Mexico 
Border and 
Gulf Coast 


West Coast 

u.s. 


90°W 


22 March 
1200 GMT 


U.S, -Canada 
Border 


U.S. -Mexico 
Border and 
Gulf Coast 


115°W 


80°W 


23 March 
0000 GMT 


U. S. -Canada 
Border 


U.S, -Mexico 
Border and 
Gulf Coast 


105°w 


75°W 


23 March 
1200 GMT 


U.S. -Canada 
Border 


U.S. -Mexico 
Border and 
Gulf Coast 


105°w 


75°W 


24 March 
0000 GMT 


U.S. -Canada 
Border 


U.S. -Mexico 
Border and 
Gulf Coast 


105°w 


65°W 
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SECTION III 



THE MONTGOMERY STREAM FUNCTION COMPUTER PROGRAM 
- o 

Assuming d0/dt , the equations of motion in 0 coordinates become 

A 
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where the isentropic stream function is defined as 



AA - c P T + <3 Z . ( 3 ) 

For the air trajectory method as derived by Danielsen (1961) and 
later outlined by Mahlman (1965 ) , the Montgomery stream function is 
given by 

NA — Cp"^ A ^ +■ R T~ C9f\ Fp ~ fg) . (2,.) 



As observed by Danielsen (1959), large errors of M result if val- 
ues of Tq and Pq are selected independently of each other from the 
plotted soundings. Consequently, for this program, the dependent re- 
lationship of Tq and Pq on an isentropic surface, as given in Poisson*s 
equation, allows M to be written 





(5) 



as shown in Mahlman and Kamm ( 19651 M is now in a convenient compu- 
tational form which can be handled readily by the computer. 

The basic program is presented in Appendix I and will be briefly 
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outlined in this section. After reading in the data, all temperatures 
are converted to degrees Kelvin. Potential temperatures are next com- 
puted for each level. For this study, a range of potential temperatures 
between 290 K and 350 K at 5° intervals is assigned for computation. 

The data i» now searched for the first 0 level above and below the 
predetermined 0 level. With these two 0 values and the predetermined 
0 value, Tq is computed. Pq is then computed from Poisson’s equation 
using and 0. Finding next the first standard levels above and below 
the arbitrary isentropic surface and calculating T by summing the weight- 
ed mean temperatures between all data points from the lower standard 
level to the arbitrary 0 level, the Montgomery stream function is com- 
puted in three segments according to equation (5). Finally, the wind 
direction and wind speed k ifc interpolated according to the procedure used 
by Mahlman and Kamm (1965). After completing the calculations for one 
particular isentropic surface at all stations in the data area, the pro- 
cedure is repeated for the next higher relevant 0 level. 

The final computer output gives for each isentropic surface the M 
value, pressure, static stability, wind direction, and wind speed for 
all reporting stations. With these five parameters, the isentropic 
trajectories can be computed and checked for accuracy, thus enabling 
the three-dimensional tracking of selected air parcels. 
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SECTION IV 



DETERMINATION OF RELEVANT ISENTROPIC SURFACES 

Possessing now the computed parameters M, Pq, wind direction and 
wind speed for all the isentropic surfaces at 5° intervals between 290 K 
and 350 K, the next step in the procedure is to determine particular 0 
levels which are representative of the isentropic surfaces present in 
the polar front jet and the subtropical jet streams. 

The observation time 23 March 12000 GMT is selected as the synoptic 
map most representative of the five observation times. A north-south 
cross section on longitude 95° W is constructed with plotted values 
of wind speed and potential temperature as depicted in Figure 2. An- 
alysis of this cross section shows that four isentropic surfaces might 
possibly be of interest — 310 K, 320 K, 330 K, and 335 K. ‘These four 
0 surfaces are chosen because of their physical location with respect 
to both jet streams and their heights above sea level, so that there 
will be a minimum of orographic and lower level interference. 

Basic maps are plotted on 1:15,000,000 scale charts for the four 
0 surfaces. Each of these isentropic surfaces is analysed for the 
Montgomery stream function using a 0,5x10^ crn^ sec - ^ interval, for pres- 
sure using a 50 mb interval, and for wind speed using a 10 m sec"-*- iso- 
tach interval. Figures 3, 4, and 5 are typical examples of these anal- 
yses. 
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Fig. 2, 23 March 1966, 1200 GMT north-south cross-section on longi- 

tude 95°W. dolid lines aro isotachs (units m sec -1 ); dashed lines are 
potential temperature (°K). 310 K, 320 K, 330 K, and 335 K are the 

isentropic surfaces chosen for trajectory analysis in this study. 
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310 K isotach analysis for 23 March 1966, 1200 GMT, Wind speed expressed in m sec 
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310 K pressure analysis for 23 March 1966, 1200 GMT. Pressure expressed in mb 



SECTION V 



COMPUTATION OF ISfiNTROPIC AIR TRAJECTORIES 



Mahlman’s outline (1965) of the Danielsen method (1961) develops 
in detail the two working equations for computation of isentropic air 
trajectories , 



The left hand term in equation (6) represents the local change in 
the stream function. The local change is graphically determined by an- 
alysing the 12 hour local changes of nM at each upper-air station or 
data point. For the four time periods of this study, &M is analysed 



trajectory path to obtain an approximation to the left hand side of 
equation (6). Figure 6 illustrates the local change of M for the period 
22 March 0000 GMT to 22 March 1200 GMT on the 330 K 0 surface. 

The second term on the right hand side of equation (6) expresses 
the change in kinetic energy of the horizontal wind. After converting 
the units of this kinetic energy term to the units of M (10? cm^ sec”^) 
and ra-arranging terms, equation (6) can be written as 




( 6 ) 



D = \4 At. 



(7) 



using an increment of 0.25x10 1 cm^ sec"^. aM is then averaged along the 



M 2 = M, * O.ooos (V 2 - V 2 ) a + 2fM 



( 8 ) 



for V expressed in m sec"^-. 
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The units of equation (7)> the distance equation, are converted to 
degrees latitude for a 6 hour trajectory, becoming 

DiO./9V£. (9) 

The objective now is to locate a value of M at the final point on 
a trajectory path such that equations (8) and (9) are simultaneously 
satisfied,, Detailed procedures for calculating isentropic air traject- 
ories can be found in Danielsen's paper (1961). However, the procedures 
used for this study will be briefly outlined for one time period. 

For the first 12 hour period, 22 0000 GMT to 22 1200 GMT, a rep- 
resentative sample of initial data points in and around the isotach 
center of the polar front jet and the subtropical jet streams is selected 
on the 22 0000 GMT M analysis for the four relevant 9 surfaces. The 
first estimate of the path of the trajectory is obtained by sketching 
both the 22 0000 GMT and the 22 1200 GMT streamlines around one of the 
initial data points. An average streamline is determined by giving 
more weight to the 22 0000 GMT streamlines for the first 6 hour period 
and to the 22 1200 GMT for the second 6 hour period. A 6 hour traject- 
ory is computed by equation (9) using the wind speed at the initial data 
point and is plotted along the average streamline using 1° increments. 

A mean wind at the 6 hour position is determined from the 22 0000 GMT 
and 22 1200 GMT isotach analyses. This mean wind is then averaged 
with the initial wind to determine an estimated mean wind over the 
first 6 hour period. Equation (9) is again employed to calculate a 
better estimate of the 6 hour trajectory position. 

Using the mean wind at the 6 hour position, the trajectory for the 
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next 6 hour period is computed by equation (9). The winds at the 12 
hour and 6 hour positions are averaged and this resulting mean wind is 
used to re-estimate the 12 hour position. 

Now, the trajectory computed by using equation (9) must be tested 
by equation (8), To, calculate AM, the trajectory determined by the 
distance equation is superimposed on the respective aM analysis and 
divided into four equal segments, Simpson's rule is utilized to calcu- 
late the 12 hour local change of M integrated along the trajectory path. 
Having aM, Vj_, V £, and M-^, M 2 is computed. 

At this stage in the procedure, an observed value of M 2 , determined 
from the first estimate to the trajectory and the 22 1200 GMT M analysis, 
and a computed value of M 2 , determined by equation (8), have been ob- 
tained, M2_obs com P arec * ^2-comp* ^ their absolute difference 

is equal to or less than 0,1x10* cm^ sec”^, the initial trajectory sat- 
isfies both equations and this trajectory represents the 12 hour air 
parcel movement on the isentropic surface. 

Generally, for the first estimated trajectory, j ^2- C omp “ ^2-obs| 
is greater than 0,1x10^ cm^ sec“^. If this is the case, the 12 hour 
position calculated by equation (9) is moved perpendicular to the mean 
flow towards the M 2 _ com p value obtained from equation (8), A new V 2 at 
this second estimated 12 hour position is determined from the 22 1200 
GMT isotach analysis. 

The initial trajectory is now adjusted throughout its entire path 
to conform to the second estimated 12 hour position. The resulting 
differences between the adjusted trajectory and the initial trajectory 
have changed the 6 hour position. Hence, a new mean wind over the first 
6 hour period and the second 6 hour period is calculated by the method 
outlined previously. 
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Following the original procedures with the new 6 hour mean winds, 
a second estimated 12 hour trajectory position on the adjusted trajec- 



tance and the energy equations simultaneously. If not, the previous 
procedures are again repeated until convergence is reached. lhus, an 
air parcel has been successfully tracked for 12 hours on an isentropic 
surface by successive approximations. 

Figure 7 illustrates the final estimated trajectories for the 
22 0000 GMT to 23 0000 GMT time period on the 330 K isentropic surface. 



tory is computed. If is less than or equal to 

0.1x107 cm^ sec“2 ? the second estimated trajectory satisfies the dis- 
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Fig. 6. 310K uM analysis for 22 March 1966, 0000 GMT to 1200 <34T. aM expressed i 
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SECTION VI 



RESULTS OF THE ISENTROPIC AIR TRAJECTORY COMPUTATIONS 

After computing a representative sample of trajectories for each 
of the four isentropic surfaces, (310 K, 320 K, 330 K, and 335 K), for 
all four time periods, the trajectories are superimposed on the respec- 
tive pressure analyses, so that Pq can be determined for the initial and 
final position, £>is then calculated for each 12 hour trajectory, where 



To illustrate the three-dimensional air parcel movement, an analy- 
sis of trajectory position relative to the surface front position is 
made. The positional relationship between a trajectory and the moving 
surface frontal system, as shown in Figure 7, appears quite consistent 
for all four 0 surfaces. Trajectory 2 on the 330 K surface, initially 
over Las Vegas, Nevada (Station 386), has an average downward vertical 
motion (<*J positive) of +16 mb/12hrs for the 22 0000 GMT to 22 1200 GMT 
period. This downward vertical motion is reasonable since the Trajec- 
tory 2 position is generally lagging the cold front position for the 
12 hour period — thus remaining in a subsidence region. 

As Trajectory 2 progresses to the northeast, its position relative 
to the frontal system is, on the average, leading the surface system. 
This leading situation is indicative of the average upward vertical 
motion of -1 mb/12hrs for the 22 1200 GMT to 23 0000 GMT period. 

This leading and lagging relative relationship between isentropic 
air trajectories and the surface frontal system was consistent for all 
time periods and isentropic surfaces. Thus, through the laborious 4 



- r tL 

CaJ — Jt 

At J d t 

*1 




( 10 ) 
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level computations, a physical understanding of air trajectory three- 
dimensional motion has been attained. In addition, a mental picture of 
this motion can be generated by relating these trajectories to the sur- 
face phenomenon. 
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SECTION VII 



TESTING THE HYPOTHESIS 

Recent studies by Danielsen (1961), Reiter and Whitney (1965), and 
Manlman (1965) have investigated cases of apparent jet stream cross-over 
at isentropic levels below 310 K. Assuming the cross-over phenomenon 
does exist at the lower levels of the atmosphere, it is reasonable to 
inquire as to the existence of a compensatory cross-over flow above 
the 310 K isentropic level. Therefore, this hypothesis will be tested 
by further analyses of trajectories associated with the respective jet 
streams within the merging region. 

The isentropic air trajectories are analysed for their positions 
relative to isotach centers of the polar front jet and the subtropical 
jet streams. The trajectories are further analysed for their positions 
relative to the region of apparent jet stream merging. As illustrated 
in Figure 3, Trajectory 1 at 22 0000 GMT is located within the 50 m 
sec - '*' isotach center of the polar front jet stream on the 335 K sur- 
face, Maintaining a cyclonic curvature with an average downward vertical 
motion of +8 mb/12hrs for the succeeding 12 hour period, an air parcel 
associated with Trajectory 1 remains in the polar jet flow and by 22 
1200 GMT is well out of the merging region. Referring again to Figure 8, 
an air parcel associated with Trajectory 2 follows a similar path as 
that of Trajectory 1, but is displaced more to the south throughout 
the 22 0000 GMT to 22 1200 GMT period. As before, this trajectory can 
be tracked through the merging region and on to the northeast. 

After analysing the trajectories associated with the polar front 
jet stream on the four 0 surfaces for the four time periods, the mean 
three-dimensional flow of air parcels in the polar front jet is deter- 
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mined to be basically cyclonic in nature — descending when moving south- 
east and ascending when moving northeast. This mean flow corroborates 
previous studies. 

For the 22 0000 GMT to 22 1200 GMT period, the apparent jet merging 
region is over New Mexico and Arizona. Trajectory 3 of Figure 8 at 
22 0000 GMT begins within the merging region over Tucson, Arizona 
(Station 274). An air parcel following this trajectory appears to 
disassociate from the main flow of the subtropical jet and proceeds to 
the northeast with an average upward vertical motion of -10 mb/12hrs 
for the first 12 hour period. By 22 1200 GMT, the air parcel associated 
with Trajectory 3 has entered the isotach center over Missouri and 
continues to move northeastward toward the Great Lakes region. This 
phenomenon of apparent disassociation of air parcels from the subtropical 
stream flow is evident on the four 0 surfaces for the trajectories with 
initial position over Tuscon at 22 0000 GMT, 

Trajectory 4 of Figure 8 is located over El Paso, Texas (Station 
270) at 22 0000 GMT. Unlike Trajectory 3, an air parcel associated with 
Trajectory 4 remains within the anticyclonically curving subtropical 
jet stream flow for the succeeding 24 hour time period. 

To substantiate the statement made in the previous paragraph con- 
cerning apparent disassociation of air parcels from the subtropical jet 
stream flow on a 0 surface other than 335 K, refer again to Figure 6. 
Trajectory 3 on the 330 K surface also begins over Tucson, Arizona at 
22 0000 GMT. This trajectory path appears to cross-over and by 22 1200 
GMT possesses the characteristic flow pattern of the polar front jet 
stream. 

By 23 March 0000 GMT, the jet stream merging region is located over 

Texas, Trajectory 2 on the 330 K isentropic surface in Figure 9, 
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initially located over Del Rio, Texas (Station 261), appears to continue 
within the subtropical jet influence and maintains an average upward 
vertical motion of -29 mb/12hrs for the 23 0000 GMT to 23 1200 GMT 
time period. 

Ihe isotach analyses of the four isentropic surfaces for the 23 
1200 GMT to 24 0000 GMT period define* the merging region of the polar 
front and the subtropical jet streams most clearly. Trajectory 3 of 
Figure 10 was selected because not only was the air parcel following 
this trajectory path associated with the subtropical jet flow, but the 
23 1200 GMT position over Midland, Texas (Station 265) was well within 
the merging region. Trajectory 3 on the 320 K and 330 K isentropic 
surfaces appears to cross-over and proceed cyclonically to the north- 
east, apparently identifying itself with the polar front jet flow at 
these levels as illustrated in Figure 10. However, trajectories on the 
335 K and 310 K surfaces located over Midland at 23 1200 GMT are dis- 
placed to the south of the 320 K and 330 K trajectory position, indicating 
a continued association with the subtropical jet stream flow. It is 
interesting to note that Trajectory 4, located initially on the anti- 
cyclonic side of the subtropical jet stream, apparently remains well 
within this flow. 

The subtropical jet is located north of its normal position for the 
23 March 1200 GMT to 24 March 0000 GMT period. Trajectory 3 of Figure 
10 is initially located on the cyclonic side of the jet and in the 
central section of the merging region. Trajectories computed to the 
south, or on the anticyclonic side of the subtropical jet stream, follow 
the classical anticyclonic path over the southern United States as 
shown in Figures 7 and 9. 
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Potential vorticity (P) of an air parcel is known to be a conserva- 
tive property for adiabatic and frictionless air motions. As a further 
test for the four trajectories with apparent cross-over flow, P is com- 
puted for the initial, 12 hour, and 24 hour isentropic trajectory posi- 
tions. If the three P values do not vary appreciably, this lends further 
credence to the validity of the trajectories indicating cross-over flow. 
Consequently, the fact that potential vorticity is conserved, as shown 
in Table 2, helps substantiate the results of the isentropic trajectory 
computations , 

v 

As shown in Figure 8, Trajectory 3 on the 335 K isentropic surface 
for the first 12 hours follows a cyclonic path and for the second 12 
hours an anticyclonic path. If the assumption of conservation of po- 
tential vorticity is valid for this assumed frictionless and adiabatic 
motion, then the large 24 hour relative vorticity change which reflects 
this curvature change must be explained by a significant increase of 
one of the parameters in the P equation. One possibility is a rapid 
increase in the Coriolis parameter, but for this trajectory, variations 
in f are small as shown in Table 2. However, Table 2 shows that the 
relative vorticity decrease is compensated by a rapid static stability 
increase. To some extent this argument helps justify the physical 
mechanisms involved in the cross-over phenomenon. 
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Fig. 8. Trajectories on the 335 K isentropic surface. 22 March 0000 GMT isotach centers, solid 
lines; 22 March 1200 GMT isotach centers; dashed lines. Hatched area shows approximate jet merging 
region. Refer to rig. 6 for units of M, Pq, and dd. © (mb/12hrs) entered above each trajectory. 
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Table 2. Approximate values of potential vorticity computed at rela- 
vant trajectory positions. Note the static stability and absolute vor- 
ticity values corresponding to respective trajectory curvatures shown 
in Figures 7 , 3 , and 10. 



0 SURFACE 





330 K 


335 K 


320 K 


1 330 K 

I 


FIGURE 


7 


8 


10 


10 


TRAJECTORY 


3 


3 


3 


3 




22 0000 GMT 


22 0000 GMT 


23 1200 GMT 


23 1200 GMT 


DATE 


22 1200 GMT 


22 1200 GMT 


24 0000 GMT 


24 0000 GMT 




23 0000 GMT 


23 0000 GMT 






STATIC 


3. 8x10" ^ 


3. 8x10" 5 


0.53x10" 5 


5. 6x10" 5 


STABILITY 


4. 5x10" ^ 


5.4xlO"5 


0.50x10"5 


5. 3x10" 5 


(deg cm sec + ^g" 4 ) 


S.OxlO" 5 


8 o 0xl0"5 






CORIOLIS 


0.73xl0" 4 


0.73x10"^ 


0.79xl0" 4 


0.79xl0" 4 


PARAMETER 


0.83xl0" 4 


0.83xl0" 4 


0.87xl0" 4 


0.83xl0" 4 


(sec"- 1 -) 


0,93xl0" 4 


0.93xl0" 4 






SHEARING 


0.30x10t4 


0.30xl0" 4 


-0.20xl0~ 4 


0.0 


VORTICITY 


o.o 


0.12xl0" 4 


-0.25xl0" 4 


-0.06xl0" 4 


(sec" 4 ) 


0.0 


-0.06xl0" 4 






CURVATURE 


0.28xl0" 4 


0.27xl0" 4 


0.27xl0" 4 


0.32xl0" 4 


VORTICITY 


0.27x10"^ 


0.3Ixl0" 4 


0.34xl0" 4 


0.30xl0" 4 


(sec"- 1 -) 


-0.28xl0" 4 


-0.22xl0" 4 






ABSOLUTE 


1.31xlO" 4 


1.30x10"^ 


0.86xl0" 4 


l.llxlO" 4 


VORTICITY 


l.lOxlO" 4 


l.lOxlO" 4 


0.96x10" 4 


1.19xl0" 4 


(sec"- 1 -) 


0 o 65xl0“ 4 


0.65xl0" 4 






POTENTIAL 


5xlCr9 


5x10" 9 


5x10"9 


6x10" 9 


VORTICITY 


5xlO" 9 


6xlO" 9 


5xlO"9 


6xlO"9 


(deg cm sec g" 1 ) 


5xlO“ 9 


5x10" 9 







40 



SECTION VIII 



CONCLUSIONS 

As indicated in figures 7, 8, and 10 and their subsequent analyses, 
cross-over of air parcels in the subtropical jet stream flow to the 
polar front jet flow apparently exists between the 310 K and 335 K 
isentropic surfaces. This conclusion is dependent upon the validity of 
the adiabatic assumption. However, ref erring again to Figures 7, 8 
and 10, it can be observed that small directional wind shears in the 
vertical exist through the 10° isentropic layer. Therefore, diabatic 
effects, even if present, would not appreciably affect the results. An 
additional implication from these analyses is the existence of a mass 
transport phenomenon in the jet stream region as observed by Reiter and 
Whitney (1965). 

The cross-over phenomenon of air parcels initially in the polar 
front jet stream is not observed in this case study for the relevant 
isentropic surfaces. However, the apparent lack of air parcel cross- 
over to the subtropical jet stream flow at these 0 levels might be ex- 
plained by the limited number of isentropic surfaces, synoptic situa- 
tions, and parameters utilized in this research. 
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APPENDIX I 



The Montgomery stream function computer program developed for this 
study (Sections A and B). 
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APPENDIX I; Section A, List of symbols used in the computer program 



SYMBOL 

DD1 

DD2 

DDTH 

F 

FI 

F2 

F3 

FF1 

FF2 

FFTH 

PP 

PP2 

PPB 

PPT 

PTH 

STAB 

T 

TB 

TBAR 

TH 

THB 

TH1 

THT 

TT 



MEANING 

wind direction at first standard level below TH1 

wind direction at first standard level above TH1 

interpolated wind direction at TH1 

Montgomery stream function 

first segment of F 

second segment of F 

third segment of F 

wind speed at first standard level below TH1 
wind speed at first standard level above TH1 
interpolated wind speed at TH1 
pressure at first standard level below TH1 
pressure at first standard level above TH1 
pressure at first level below TH1 
pressure at first level above TH1 
pressure at TH1 



stability parameter = - 



ae_ 

9P 



temperature 

T at first level below TH1 

average T of layer between PP and PP2 

potential temperature 

TH at first level below TH1 

assigned TH value 

TH at first level' above TH1 

T at first level above TH1 
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3IMB0L 



MEANING 



TTH 

ZTH 



T at 1H1 

geopotential height of TH1 
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onoooo 



APPENDIX I; Section B. The Montgomery etream function computer program. 



MONTGOMERY STREAM FUNCTION COMPUTER PROGRAM — 1 JUNE I960 

RADIOSONDE DATA FOR 22 MARCH 0000 GMT TO 24 MARCH 0000 GMT 
1966. FIVE OBSERVATION TIMES -- RESPECTIVE MATRIX A,B,C,D,E. 



DIMENSION A< 1190* 5 ) , 
1 G ( 1484, 5 ) » THI1484) 



500 

600 



READ 

READ 

READ 

READ 

READ 



(5,500) 
(5,500) 

( 5,500) 
(5,500) 

(5.500) 

FORMAT (10F8.0) 
WRITE (6,600) 
FORMAT (IH1) 

N = 1 

.EQ. 

• EQ. 

.EQ. 

.EQ. 

.EQ. 

.EQ. 



m ihohi , 

( ( A( I , J) , 
( ( B ( I , J ) , 
( ( C ( I , J ) , 
( ( D ( I , J ) , 
( ( E ( I , J ) , 



IF 

IF 

IF 

IF 

IF 

IF 

M 

GO 



(N 
( N 
(N 
(N 
(N 

(N ... 

= 1190 
TO 7 



1) 

2) 

3) 

4) 

5) 

6 ) 



GO 

GO 

GO 

GO 

GO 

GO 



M = 1478 
GO TO 7 
M = 1270 

GO TO 7 
M * 1290 
GO TO 7 
M = 1484 



B( 1478,5) , 

T ( 1484) 

J = 1,5). 

J = 1, 5), 
J * 1, 5) , 

J = 1,5), 

J = 1,5), 



0(1290,5), E(1484,5>, C(1270,5) 



TO 

TO 

TO 

TO 

TO 

TO 



7 


DO 


10 


i = : 


L, 




M 








DO 


10- 


j = 


1, 




5 








IF 


(N 


.EQ. 


1 


) 


GO 


TO 


20 




IF 


( N 


.EQ. 


2 


) 


GO 


TO 


21 




IF 


(N 


.EQ. 


3 


) 


GO 


TO 


22 




IF 


(N 


.EQ. 


4 


) 


GO 


TO 


23 




IF 


(N 


.EQ. 


5 


) 


GO 


TO 


24 


20 


G ( I 


,J) 


= A 1 


( I 


f 


J) 








GO 


TO 


10 












21 


G( I 


, J ) 


= B 


( I 


t 


J) 








GO 


TO 


10 












22 


G( I 


,J) 


= C 


( I 


f 


J) 








GO 


TO 


10 













I = 1, 1190) 

I = 1, 1478) 
I = 1, 1270) 

I = 1, 1290) 

I = 1,1484) 



2 

3 

4 

5 

6 

17 
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non non non 



23 G( I.J) * D< I , J) 

GO TO 10 

2 A G( I * J) = E ( I » J ) 

10 CONTINUE 

25 DO 11 I = 1, M 

11 WRITE (6.610) I, ( G( I . J ) * J = 1, 5) 

610 FORMAT (20X, 110, 5X , 5F10.1) 

CONVERT TEMPERATURES AND COMPUTE POTENTIAL TEMPERATURES 

WRITE (6,612) 

612 FORMAT ( lHl ) 

DO 12 I = 1, M 

IF ( G( 1,3) .E0. 777.0) GO TO 8 
T( I) = 6 ( 1 , 2 ) + 273.2 

TH ( I ) = T ( I ) * ((1000.0 / G( 1 , 1 ) ) ** 0.2857) 

GO TO 18 
e T ( I ) = 777.0 
TH( I ) = 777.0 

18 WRITE (6 ,611 ) I , T( I ) , TH( I ) 

611 FORMAT (5X, 110, 2F10.1) 

12 CONTINUE 

CHECK FOR DI SCREPENCIES IN POTENTIAL TEMPERATURES ABOVE 800 MB 
WRITE (6,620) 

620 FORMAT (1H1, T10, • D I SCREP ENC I ES IN THETA VALUES', ///) 

DO 13 I = 1, M 

IF (TH(I) .EQ. 777.0) GO TO 13 
IF (G( I , 1) - 800.0) 40, 40, 13 

40 IF ( TH( 1+ 1) - TH ( I ) ) 41, 13, 13 

41 WRITE (6,630) I, TH(I), THlI+1) 

630 FORMAT <$X, 110, 2F10.1, /// ) 

13 CONTINUE 

ASSIGN POTENTIAL TEMPERATURES FOR COMPUTATION 

TH 1 = 290.0 
GO TO 50 

14 TH 1 = TH1 + 5.0 

IF ( TH1 .GT. 350.0) GO TO 16 
50 WRITE (6, 640) 

640 FORMAT ( lHl ) 

FF TH = 0.0 
DOTH = 0.0 
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SEARCH DATA FOR VALUES ABOVE AND BELOW LEVEL BEING COMPUTED 



DO 15 I = 1, M 

IF ( TH( I I .EQ. 777.0) GO TO 15 
IF ( TH ( I i - TH1 ) 51, 52, 53 

51 GO TO 15 

52 TTH = T(I) 

GO TO 54 

53 THT = TH ( I ) 

TT = T ( I ) 

PPT = G( 1,1) 

IF (TH(I-l) .EO. 777.0) GO TO 15 
IF ( TH( I - 1 ) - TH1 ) 55, 15, 15 
55 TB = T ( I — 1 ) 

PPB = G< 1-1,1) 

THB = THU-1) 

COMPUTE PRESSURE AND TEMPERATURE AT LEVEL OF INTEREST 

TTH = TB + (TT - TB) * UTH1 - THB) / (THT - THB)) 

54 PTH = 1000.0 * ((TTH / TH1 ) ** 3.5001) 

FIND NEAREST STANDARD LEVEL 

IPP = PTH / 50.0 
AAA * IPP * 50 
L = I 

60 IF (G(L,1) .EO. AAA) GO TO 61 
L = L + 1 

GO TO 60 

61 PP = G ( L , 1 ) ♦ 50.0 
J = L 

62 IF (G(J,1) .EO. PP) GO TO 63 
J = J - 1 

GO TO 62 

63 PP = G ( J , 1 ) 

HGT = G ( J, 3 ) 

FF1 = G( J, 5 ) 

DD1 = G(J,4) 

AA = 0.0 

70 IF (PTH - G ( J+l , 1 ) ) 71, 72, 72 

71 AA = AA + ( T ( J I + T ( J+l ) ) * (G(J,1) - G(J+1,1)) 

J = J + 1 

GO TO 70 



U8 



72 AA = AA ♦ ( T ( J ) ♦ TTH) * (G(J,1) - PTH) 

74 TBAR = (1.0 / (2.0 * ( PP - PTH))) * AA 
PP2 - PP - 50.0 

6»* IF ( G ( J « l ) .E0. PP2) GO TO 65 

J = J ♦ 1 
GO TO 64 
65 PP2 = G( J, 1 ) 

FF2 = G ( J » 5 ) 

002 = G(J,4) 

C 

C COMPUTE MONTGOMERY STREAM FUNCTION IN THREE SEGMENTS 

C 

Z = 1.E6 

FI = 10.046 * TTH * Z 
F2 = 9.8 * HGT * 1.E4 

F3 = 2.8704 * Z * TBAR * (ALOG(PP) - ALOG( PTH) ) 

C 

C INTERPOLATE wINO DIRECTION AND WIND SPEEO 

C 

IF (FF1 .IT. 8B8.0) GO TO 80 
FFTH = 888.0 
DOTH = 888.0 

8C IF (FF2 .LT. 888.0) GO TO 81 
FFTH = 888.0 
DOTH = 888.0 
GO TO 88 

81 FFTH = FF1 * IFF 2 - FF1 ) * ((PP - PTH) / 50.0) 

DO 3 = DD2 - DO 1 

IF ( DD3 - 180.0) 82, 83, 83 

83 DD3 = DD3 - 360.0 
GO TO 85 

82 IF ( 0D3+ 180.0) 84, 85, 85 

84 DD3 = DD3 ♦ 360.0 

85 DOTH = DD1 ♦ DD3 * ( ( PP - PTH) / 50.0) 

IF (DOTH) 86, 88. 87 

86 DOTH = DOTH ♦ 360.0 
GO TO 88 

87 IF (DOTH - 360.0) 89, 89, 89 
89 DOTH = DDTH - 360.0 

88 CONTINUE 
C 

C COMPUTE STABILITY FACTOR 

C 

STAB = ( THT - th«) / (PPB - PPT) 
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COMPUTE MONTGOMERY STREAM FUNCTION 

F = FI + F2 ♦ F3 

COMPUTE HEIGHT OF THETA SURFACE 

ZTH = HGT - (2.8704 * TBAR / 9.8) * (ALOG(PTH) 

WRITE (6,650) I, STAB, FI, F2, F3, F 
650 FORMAT (5X, 110, 5X, F6.3, 5X, 3P4E15.5) 

WRITE (6,660) TH1 , PTH, TTH , ZTH, DOTH, FFTH 
660 FORMAT (5X, 6F10.1, ///) 

15 CONTINUE 
GO TO 14 

16 N = N + 1 
GO TO 1 

17 STOP 
END 



ALOG(PP) ) 
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